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THE  INVERSE  SOURCE  PROBLEM  FOR  MAXWELL’S  EQUATIONS 

R.  ALBANESE*  AND  P.B.  MONK* 

Abstract.  The  inverse  source  problem  for  Maxwell’s  equations  is  considered.  We  show  that  the  problem  of  finding  a 
volume  current  density  from  surface  measurements  does  not  have  a  unique  solution,  and  characterize  the  non-uniqueness.  We 
also  show  that  if  further  a  priori  information  is  available,  the  inverse  source  problem  may  have  a  unique  solution  (in  particular 
for  surface  currents  or  dipole  sources). 

1.  Introduction.  The  goal  of  this  paper  is  to  investigate  the  inverse  source  problem  for  the  general 
Maxwell  system.  This  problem  arises  in  medical  applications  where  measurements  are  taken  of  the  electric 
and  magnetic  surface  currents  on  the  surface  of  the  human  head,  and  it  is  desired  to  infer  from  these  currents 
the  source  currents  in  the  brain  that  produced  the  measured  fields.  It  is  hoped  that  such  measurements  could 
be  used  to  diagnose  abnormalities  in  the  brain  and  also  to  allow  the  control  of  prosthetic  limbs. 

From  the  point  of  view  of  mathematical  modeling,  we  shall  first  make  the  simplifying  assumption  that 
the  measurements  can  be  made  on  a  surface  containing  the  entire  head.  Later  we  shall  discuss  how  the 
theory  developed  in  this  ideal  case  might  be  extended  to  the  more  realistic  problem  where  measurements  are 
made  on  a  portion  of  the  skull.  Thus  we  assume  that  there  is  a  bounded  smooth  domain  Q  C  R3  (connected 
with  connected  complement)  such  that  the  known  electromagnetic  parameters  e,  a  and  p  (all  assumed  real) 
have  the  following  properties: 

1.  In  R3  \  Q  (the  air)  the  parameters  have  the  values  e  =  eo  >  0  and  p  =  po  >  0  where  eo  and  Ho  are 
the  permitivitty  and  permeability  of  free  space  respectively.  In  addition  the  conductivity  cr  =  0. 

2.  In  0,  (the  head)  there  are  constants  tmin  and  emax  such  that 

0  ^  £ min  ^  e(3?)  5^  £max  fo7  all  X  £  fl. 

In  addition  there  are  a  constant  omax  such  that 

0  <  a{x)  <  amax,  for  all  x  £  0. 

In  view  of  the  biological  application  it  is  also  reasonable  to  assume  that  p  =  fio  in  all  O  and  we  shall 
make  this  assumption  in  the  rest  of  the  paper. 

3.  The  coefficients  e  and  a  are  assumed  to  be  piecewise  C1  functions  and  any  surfaces  of  discontinuity 
are  assumed  to  be  smooth. 

The  latter  assumption  is  made  so  that  we  can  apply  the  unique  continuation  results  of  [21]  (in  [19]  it 
is  commented  that  the  continuation  property  extends  to  piecewise  C1  functions  under  even  more  general 
conditions  on  the  interface  between  regions  where  the  coefficients  are  smooth). 

Given  a  current  density  J  supported  in  SI,  the  resulting  electromagnetic  field  described  by  the  electric 
field  E  and  magnetic  field  H  satisfies  the  following  Maxwell  system  where  w  >  0  is  the  frequency  of  the 
time  harmonic  field: 

(1.1)  -\ueE  +  aE  -  V  x  H  —  -J  in  R3, 

(1.2)  — iwpff  +  V  x  E  =  0  in  M3. 

In  order  to  uniquely  determine  the  fields,  we  must  also  impose  the  Silver-Miiller  radiation  condition  at 
infinity: 

(1.3)  lim  ( H  xx  —  \x\E)  =  0 

|ac  | — >00 

uniformly  in  x  =  x/\x\.  Under  suitable  conditions  on  the  smoothness  of  the  source  J  (for  example  J  £ 
L2(R3)  having  compact  support),  the  forward  problem  of  finding  E  and  H  for  a  known  current  source  J  is 
well  posed  (see  for  example  [19]). 
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Having  defined  the  above  notation,  we  can  now  discuss  the  inverse  source  problem.  Let  T  denote  the 
surface  of  il.  Let  n  denote  the  unit  outward  normal  to  V  =  (Kl.  For  the  inverse  source  problem  we  assume 
that  e,  a  and  /r  are  known  and  in  addition  that  the  surface  electric  current  n  x  E  and  the  surface  magnetic 
current  n  x  H  are  known  on  the  entire  surface  T.  From  this  data  we  would  like  to  find  the  corresponding 
unknown  source  current  density  J .  As  we  shall  see  this  problem  is  severely  ill-posed. 

The  inverse  source  problem  has  been  exhaustively  studied  in  the  literature  both  from  the  point  of  view  of 
applied  biomedical  engineering  and  also  as  a  mathematical  problem  (see  for  example  [3,  22,  15,  7,  14,  9,  6,  12] 
where  we  have  emphasized  mathematical  references  of  relevance  to  this  paper).  It  is  well-known  that  volume 
sources  cannot  be  uniquely  determined  from  surface  measurements  due  to  the  existence  of  non-radiating 
sources  [3,  18,  17,  8].  One  goal  of  this  paper  is  to  extend  the  results  of  [3]  and  characterize  the  non-radiating 
sources  when  the  background  media  is  non-constant  as  is  found  in  the  head.  This  is  done  using  weak  solutions 
and  variational  methods. 

With  the  exception  of  the  work  of  Aminari,  Bao  and  Flemming  [1]  and  He  and  Romanov  [13]  most 
theoretical  work  on  biological  applications  has  focused  on  the  static  (or  quasi-static)  case  since  it  is  clear 
that  the  frequency  of  the  radiation  is  very  low  (perhaps  100Hz).  Even  in  [1]  the  numerical  scheme  suggested 
is  based  on  a  static  model.  However,  since  the  goal  of  the  inverse  source  problem  is  to  monitor  dynamic 
neuronal  events  (an  action  potential  has  a  rise  time  on  the  order  of  0.5  milliseconds  [23])  it  may  be  that  the 
displacement  current  is  not  negligible.  This  has  already  been  pointed  out  in  studies  of  source  problems  related 
to  monitoring  neurons  in  the  arm  [16],  In  particular  there  are  two  important  dimensionless  parameters  that 
must  be  considered  when  assessing  whether  the  displacement  current  may  be  neglected  [11]  namely 

o  =  s/ufjMjL  and  (3  =  ,/juouL 

where  L  is  a  representative  length  for  the  conductor.  Of  course  neglecting  the  displacement  current  results 
in  an  eddy  current  model  -  the  first  step  to  a  static  model 

If  we  choose  u>  =  100 Hz  then,  using  the  Gabriel  database  [10],  we  find  that  for  grey  matter  the  relative 
permitivity  is  cr  =  3.9  x  10°  and  a  =  O.OSOSm-1.  Taking  Co  =  8.85  x  10~12Fm_1,  fio  =  47T  x  10~7Hin-1  and 
taking  as  a  representative  length  L  =  0.01m  (roughly  the  radius  of  the  head)  we  obtain 

a  =  3.3  x  1(T5  and  (3  =  6.6  x  1CT6. 

While  both  these  terms  are  small  compared  to  unity,  they  are  of  the  same  order  of  magnitude  and  thus  we 
cannot  a  priori  neglect  the  displacement  current  (this  is  also  the  conclusion  of  the  computational  studies 
in  [16]).  It  therefore  seems  prudent  to  consider  the  full  Maxwell  system.  The  remainder  of  the  paper 
is  devoted  to  uniqueness  questions  for  the  inverse  source  problem  for  the  Maxwell  system.  Our  work  is 
strongly  motivated  by  [12]  in  which  uniqueness  questions  are  examined  for  the  corresponding  static  problem. 

The  plan  of  this  paper  is  as  follows.  In  the  next  section  (Section  2)  we  address  the  question  of  uniqueness 
for  the  inverse  problem  assuming  that  the  unknown  sources  are  distributed  throughout  the  domain  fi.  We 
show  that  such  sources  cannot  be  determined  from  the  given  data  (without  further  a  priori  information)  and 
characterize  the  non-uniqueness.  This  tells  us  the  limits  on  the  inverse  source  problem  in  the  sense  that, 
even  with  perfect  data,  the  volume  current  density  cannot  be  reconstructed.  If  the  source  currents  are  a 
priori  known  to  have  smaller  support  than  the  entire  domain  f2.  it  is  sometimes  possible  to  prove  that  the 
data  uniquely  determines  the  source.  We  consider  two  cases  found  in  applications.  In  Section  3  we  suppose 
that  it  is  known  a  priori  that  the  source  current  is  given  by  a  delta  function  that  is  supported  on  the  surface 
of  a  smooth  subdomain  in  f l.  In  the  case  of  a  surface  current  we  show  that  the  current  density  is  uniquely 
determined  by  the  data.  In  the  second  case,  in  Section  4,  we  consider  the  case  when  the  source  current  is 
due  to  finitely  many  dipole  sources.  In  this  case  we  prove  that  the  data  uniquely  determines  the  number, 
position  and  strength  of  the  dipoles.  Finally,  in  Section  5,  we  show  how  the  variational  characterization  can 
be  extended  to  allow  for  measurements  only  on  a  portion  of  the  boundary  F  and  draw  some  conclusions  in 
Section  6. 

The  uniqueness  theory  in  this  paper  does  not  address  the  stability  of  the  inverse  problem  to  data 
perturbations.  We  also  do  not  address  the  question  of  what  algorithm  could  be  used  in  each  of  the  three 
cases.  However  the  variational  characterization  of  the  current  source  that  we  shall  derive  can  be  used  as 
part  of  an  inversion  scheme. 


2 


In  general  vector  quantities  are  denoted  by  boldface  type  (e.g.  E  =  (Ei,E2,  ES)T  as  already  used  above). 
Much  of  our  theory  makes  use  of  the  standard  energy  space  for  solutions  of  Maxwell’s  equations 

H (curl;  ft)  =  {u  £  L2(ft)  |  V  x  u  £  L2(ft)}  . 

2.  Non  uniqueness  of  volume  currents.  Our  goal  in  this  section  is  to  derive  a  variational  equation 
relating  the  unknown  source  J  to  the  data  on  T.  This  derivation  is  motivated  by  the  work  of  Ammari,  Bao 
and  Fleming  [1]  who  derived  a  similar  expression  in  their  study  of  methods  for  identifying  point  sources.  We 
then  use  this  variational  equation  to  study  the  uniqueness  question  when  J  is  a  distributed  source  in  ft. 

Suppose  J  £  L2(ft)  has  support  in  ft  and  that  e,  a  £  L°°(K3)  satisfy  the  assumptions  noted  in  the 
previous  section.  Then  we  can  reformulate  (1.1) — (1.3)  on  ft  alone  using  a  suitable  Dirichlet-to-Neumann 
(DtN)  map  T  on  F.  In  order  to  define  the  map  T  we  need  the  trace  space  iI-1/2(Div;  T)  defined  by 

Ff-1/2(Div;r)  =  {u  €  H~1/2(T )  |  u  ■  n  =  0  on  T  and  Vr  ■  u  £  ff"1/2(r)} 

where  Vr-  is  the  surface  divergence.  Note  that  since  we  assumed  that  the  surfaces  in  this  problem  are 
smooth,  classical  definitions  of  the  surface  gradient,  curl  and  divergence  may  be  used  [4], 

Using  this  space,  the  map  T  can  be  defined  as  follows.  For  A  £  //~1//2(Div;  T)  we  define  T X  =  n  x  v 
where  ( u,v )  £  Hioc( curl;R3  \  ft)2  satisfy 

— iuieou  —  V  x  v  =  0  in  R3  \  ft, 

—iujfiov  +  V  x  u  =  0  in  R3  \  ft, 
n  x  u  =  A  on  T, 


together  with  the  Silver-Miiller  radiation  condition  (1.3).  It  follows  from  the  solvability  of  the  exterior 
scattering  problem  in  //i0C(curl;  K3  \  ft)  that  T  :  flr_1/2(Div;  T)  — >  /f~1/2(Div;  T).  For  an  introduction  to 
vector  trace  spaces  and  the  properties  of  T  see  [20,  19,  2]. 

Later,  we  shall  also  need  to  use  the  adjoint  operator  T*  :  J7-1/2(Curl;  F)  — >  jff-1/2(Curl;  T)  where 

H”1/ 2(Curl;  T)  =  {«  €  |  u  ■  n  =  0  on  r  and  Vr  x  u  £  tf-1/2^)} 

and  where  Vr  x  u  is  the  surface  curl  of  it.  Of  course  if_1/2(Div;  T)  and  Lf-1/2(Curl;  F)  are  dual  spaces  so 
that  T  and  T*  are  related  by 


J  TX  rjdA  =  J  X  T'rjdA  VA  €  /r1/2(Div; 7)  and  rj  £  IE 


2  (Curl;  T). 


The  map  T*  is  also  related  to  a  boundary  value  problem.  Given  rj  £  H  1/2(Curl;r)  define  4>  and  i/>  to 
satisfy  the  exterior  scattering  problem 

—iu>€(f)  —  V  x  i/>  =  0  in  K3  \  ft, 

—itvfiip  +  V  x  <fi  =  0  in  R3  \  ft, 

(n  x  cj>)  x  n  =  17  on  T, 
lim  (i/>  x  x  +  |as|<^)  =  0. 

|x[— »oc 

This  is  just  an  exterior  problem  for  Maxwell’s  equations  (recognizing  that  the  boundary  condition  on  T  may 
be  rewritten  as  <j>  x  n  =  17  x  n)  where  the  radiation  condition  has  the  opposite  sign  to  that  in  (1.3).  Standard 
techniques  show  that  this  problem  has  a  unique  solution  for  any  r\  £  FI_1/2(Curl;  F)  (see  for  example  [19]). 

Using  integration  by  parts  and  the  equations  for  the  fields  (it,  v)  and  (</>,  ip)  in  a  ball  Br  of  radius  R 
containing  ft  we  have  (recalling  that  n  is  the  outward  normal  to  F  and  hence  inwards  to  Br)  that 


■  /  T(n  x  u)  ■  rjdA  =  /  _  V  x  v  ■  cp  —  v  ■  V  x  <pdV  —  /  nxv-cpdA 

Jr  Jbr\U  J  dB  a 


3 


'b, ,\a 


—iuieu  ■  <p  +  V  x  u  ■  ip  dV  — 


-i u!(.u  ■  4>  +  u  ■  V  x  ip  dV 


'B„\n 

n  x  u  ■  ip  dA  + 


i 


/  : 
JdBn 


L 

l 


n  x  v  ■  4>dA 


n  x  u  ■  ip  dA  +  /  n  x  u  ■  ip  —  n  x  v  ■  cp  dA 
JdBn 


nxu-ip-nxvcp  dA. 


But,  we  may  write 

L 


/  < 

JdBn 


n  x  u  ■  —  n  x  v  ■  <p d A  =  /  (v  x  n  —  u)  ■  <p  +  u  ■  (ip  x  n  +  (p)  dA. 


Using  the  fact  that  jit|  =  0(l/|a:|),  \<p\  =  0[lj\x\),  |i;  x  n  —  u|  =  0(l/|a:|2)  and  \ip  x  n  +  <p\  =  0(l/|a:|2) 
for  [as |  =  R  large  (see  [4]),  we  see  that  the  integral  on  dBn  vanishes  in  the  limit  R  — >  oo  and  we  have  shown 
that  T*r]  —  (n  x  ip)  x  n  on  P. 

Using  the  DtN  map  T,  we  can  reformulate  (1.1 )-( 1 .3)  as  the  problem  of  finding  ( E,H )  £  //(curl;  Cl)  x 
//(curl;  il)  such  that 


(2.1) 

(2.2) 

(2.3) 


—  'uveE  +  aE  -  V  x  H  =  —J  in  fi, 
—iui/i.H  +  V  x  E  =  0  in  fi, 

H  x  n  =  T(E  x  n)  on  T. 


It,  is  now  convenient  to  derive  a  vector  wave  equation  for  E  alone  by  using  (2.2)  to  write  H  =  l/(i uifij'V  x  E 
and  eliminating  H  from  (2.1)  to  obtain 


(2.4) 


V  x  n  1 V  x  E  —  (u:2t  +  i<7o>)jE  =  iujJ  in  Cl. 


Multiplying  by  the  complex  conjugate  of  a  smooth  test  vector  function  $  and  integrating  by  parts  (using  the 
expression  for  E  just  derived)  we  obtain  the  variational  equation 


(2.5) 


//■- 


Vx£-Vx(-  (u2e  +  ia)E  ■  £  dV  +  /  iunx  H  ■  £dA=  /  iujJ  ■  £  dV. 


Jr* 


=  / 


Replacing  the  magnetic  surface  current  using  the  DtN  map  we  obtain  the  problem  of  finding  E  6  //(curl;  Cl) 
such  that 


(2.6) 


f  /u  1 V  x  E  ■  V  x  ^  —  (u;2t  +  iau>)E  ■  £  dV  +  f  iu>T(n  x  E)  ■  £  dA  =  f  iujJ  ■  £  dV 
Jii  Jr  J  n 


for  all  ^  €  H(cm\;Cl).  This  variational  formulation  can  then  be  used  to  prove  that,  for  a  given  J  £  L2(Cl), 
there  exists  a  unique  electromagnetic  field  ( E,H )  satisfying  (2.6)  [19]. 

We  now  return  to  the  inverse  problem  and  derive  an  identity  relating  the  source  current  J  to  the  data. 
The  approach  is  closely  connected  to  the  characterization  of  non-radiating  sources  in  [17]  but  we  now  allow  a 
more  general  scattering  problem  for  which  the  Maxwell  operator  is  not  self  adjoint.  Using  (2.5)  and  recalling 
that  £  is  assumed  to  be  a  smooth  vector  function,  we  integrate  by  parts  one  more  time  to  obtain  the  indentity 

iu!  f  J  ■  £  dV  =  j  E  ■  (V  x  /i-1  V  x  £  —  (uj2e  —  icrw)£)  dV 

J  n  Jn 

+  J  n  x  E  ■  fi~x  V  x  £  +  iujn  x  H  ■  £dA. 

Now  we  choose  £  £  //(curl,  Cl)  to  be  such  that 

(2.7)  /  /x_1  V  x  £  •  V  x  tp  —  (c A  -  \auj)Z  ■  ip  dV  =  0 

Jn 
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for  all  smooth  vector  functions  ip  of  compact  support  in  fi.  Note  in  particular  that  this  implies,  in  the  sense 
of  distributions,  that  £  is  just  a  weak  solution  of  the  adjoint  Maxwell  system 

V  x  /r-1V  x  £  —  ( u)2e  —  icru))£  =  0  in  fi. 

The  fact  that  £  £  L2(Q)  then  implies  that  Vx/x_1V x£  £  L2(fi)  and  hence  that  nx/r_1V x£  £  7/-1/2(Div;  T) 
is  well  defined.  Making  this  choice  of  £  we  see  that  (2.7)  becomes 

(2.8)  i  ui  J  ■  £dV  =  /  (n  x  E  ■  /x-1  V  x^  +  i un  x  H  ■  £)  dA. 

J  n  Jr 

for  all  £  £  H( curl;  fi)  satisfying  (2.7). 

Suppose  that  we  have  an  exact  knowledge  of  n  x  E  and  n  x  H  on  T  (and  of  e  and  a  in  f2).  Then 
equation  (2.8)  provides  a  link  between  the  known  data  and  the  unknown  source  J.  Let  us  define  the  set  of 
solutions  of  the  adjoint  problem  to  be 


Tt(Cl)  =  {u  £  H( curl;  Q.)  |  u  satisfies  (2.7)}  . 

Then  we  define  H{Q)  to  be  the  closure  of  7i(0)  in  the  L2(fi)  norm.  Using  H(fl)  we  can  write  the  L2- 
orthogonal  decomposition 

L2(Cl)  =  //(fi)©^^)-1. 

We  note  that  H( fi)2-  is  an  infinite  dimensional  subspace  of  L2(Cl)  as  the  following  lemma  shows. 

Lemma  2.1.  Suppose  x  €  (Co°(^))3-  Then  if 

i fi  =  V  x  x  x  ~  (w2«  —  iucr)x 

we  have  <p  £  /f(0)-L. 

Remark:  An  interesting  question  is  whether  the  above  lemma  characterizes  all  of  i/(n)-L. 

Proof.  Suppose  u  £  Ti(Q.)  then,  integrating  by  parts, 


/ 

Jn 


u  ■  <j)d.V  =  /  u  ■  (V  x  /x_1  V  x  x  -  (w2£  —  i<TW)x)  dV 
Jn 

~lV  xu-V  x  x  ~  +  i  ouj)u  ■  xdV  =  0, 


=  J  P~ 


using  (2.7).  By  a  density  argument  this  holds  for  all  functions  in  □ 

Recall  that  J  €  L2(fi).  Prom  Lemma  2.1  we  see  that  H(Q)  is  a  proper  subspace  of  L2(fl)  and  using 
(2.8)  we  see  that  only  the  component  of  J  in  H(Q)  can  be  determined  from  the  data  by  using  equation  (2.8). 
It  is  possible  that  some  other  equation  could  be  derived  to  determine  the  component  of  J  in  H{i I)2-  from 
the  data,  but  the  following  theorem  rules  this  out. 

THEOREM  2.2.  Suppose  J  £  L2(fi).  Then  J  =  J h  +  J h1-  where  J H  G  H{Sl)  and  J H±  £  The 

inverse  source  problem  uniquely  determines  J h  via  equation  (2.8).  The  component  J h ±  does  not  produce 
surface  currents  on  T  and  hence  cannot  be  identified. 

Remark:  The  functions  in  are  termed  non  radiating  sources  [3].  Computing  the  component  of 

the  source  in  H(fl)  corresponds  to  computing  a  “minumum  energy”  solution  [18]. 

Proof.  Given  the  discussion  preceding  this  theorem,  the  only  part  of  the  statement  of  the  theorem  that 
needs  to  be  proved  is  that  J H±  produces  no  surface  currents  on  T.  To  simplify  notation,  suppose  now  that 
J  £  H{Q)l.  Then  from  (2.8)  we  have  that 


L 


o 


n  x  E  ■  pT1  V  x  \  +  i  am.  x  H  ■  £  dA 
for  all  £  £  7f(Q,).  Using  the  Dirichlet  to  Neumann  map  T 

n  x  E  •  yu_1  V  x  £  +  i cjT(n  x  E)  ■  £  dA  =  0 


L 
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so  using  the  adjoint  operator  T*  :  H  1/,2(Curl;r)  — +  H  1/2(Curl;r), 

(2.9)  J  nx  E  (/F^V  x  £  -  i  wT*(£t))  <M  =  0 

where  £r  =  (n  x  £)  x  n  on  F.  We  now  choose  £  £  i/(curl;f2)  to  satisfy  (2.7)  together  with  the  boundary 
condition 


p_1(V  x  £)r  +  it uT*(£t)  =  f]  on  F 


where  rj  £  H  1/2(Curl;r).  More  precisely,  £  €  H (curl;  Cl)  satisfies  the  variational  problem 

/  /u_1  V  x  £  ■  V  x  <f>  —  (ui2e  —  io><7)£  ■  4>dV  +  j  (rj  —  iwT*(£T))  •  <fi  x  n  dA  =  0 
un  Jr 

for  all  4>  G  i/(curl;  f2).  We  have  already  shown  at  the  start  of  this  section  that  T*  is  characterized  by  solving 
an  exterior  boundary  value  problem  for  Maxwell’s  equations  using  the  radiation  condition  with  a  sign  change. 
Thus  the  above  variational  problem  is  equivalent  to  solving  the  following  strong  form  of  the  boundary  value 
problem  of  finding  £  £  7/[0C(curl;  R3)  such  that 


(2.10) 

(2.11) 

(2.12) 

(2.13) 


V  x  fi  1 V  x  £  -  (tu2e  -  icnu)£  =  0  in  R3  \  I\ 
I£J T  =  on  r, 

x  £1t  =  -rj  on  T, 
x  £)  x  x  +  itu|x|£  — »  0  as  \x\  — >  oo, 


where,  for  a  smooth  function  u  the  jump  [ii]t  on  T  is  given,  for  x  £  T,  by 

!ti]r(x)  =  lim  ((n(x)  x  u(x  -  hn(x)))  x  n(x)  -  (n(x)  x  u(x  +  hn(x)))  x  n(x))) 

h—*0 

where  we  recall  that  n  is  the  unit  outward  normal  to  Cl. 

Standard  analysis  of  uniqueness  and  existence  for  the  Maxwell  system  can  now  be  applied  to  the  above 
problem  to  show  the  existence  of  a  unique  solution  £  to  (2.10)-(2.13)  for  any  r?  £  //_1/2(Curl; Cl).  Hence 
(2.9)  shows  that 

J^nxErjdA  =  Q  for  all  r?  €  H_1/2(Curl;  Cl) 

so  n  x  E  =  0  (and  hence  H  x  n  =  T{E  x  n)  ~  0  on  T).  Thus  a  current  in  /f(ft)x  produces  no  measurable 
signal.  □ 

The  preceding  result  shows  that  only  a  certain  component  of  a  volume  current  J  can  be  determined 
directly  without  further  a.  priori  information.  Furthermore  since  the  measurable  component  is  a  solution  of 
the  homogeneous  adjoint  Maxwell  solution,  it’s  support  will  be  all  of  Cl  (unless  there  is  no  measurable  field!). 
Hence  the  support  of  the  source  also  cannot  be  determined  without  further  a  priori  information. 

3.  Uniqueness  of  surface  currents.  Next  we  consider  the  inverse  source  problem  in  a  case  where  extra 
a  priori  information  is  assumed  to  be  available.  In  particular,  we  assume  that  J  is  a  surface  current  supported 
on  a  known  smooth  surface  £  enclosing  a  domain  B  located  interior  to  Cl.  Note  that  this  assumption  is  often 
made  in  the  MEG  literature  (see  for  example  [22])  and  £  corresponds  to  the  surface  of  the  brain. 

As  we  shall  see  we  need  to  additionally  assume  that  the  surface  current  J  is  tangent  to  £  and  more 
precisely  that  J  £  //_1/2(Div;  £).  In  this  case  it  is  no  longer  true  that  H  £  H(cur\;Cl).  Instead  the 
analogue  of  equations  (2.1)-(2.3)  is  to  seek  E  £  H(curl;fi)  and  H  £  L2(Cl)  such  that  H\b  £  if  (curl;  B), 
H\n \g  €  H( curl;  Cl\B)  and 


(3.1) 

(3.2) 

(3.3) 

(3.4) 


-iujcE  +  aE  —  V  x  H  =  0  m  Cl  \  B  and  B, 
—  ic jfiH  +  V  x  E  =  0  in  fi, 

H  x  n  =  T{E  x  n)  on  T, 
[H  x  n]  =  —  J  on  £, 
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where  n  is  the  outward  normal  to  E  and  the  jump  [H  x  n], is  defined  by  [H  x  n]  =  (H |nyg  -  H\b)\s  x  n. 

The  above  equations  can  be  reduced  to  a  variational  problem  on  Si  using  the  Dirichlet-to-Neumann  map 
T.  First  we  can  eliminate  H  on  f l  \  B  and  B  respectively.  Then  multiplying  the  resulting  electric  field 
equations  by  a  smooth  test  vector  and  integrating  over  the  respective  domains,  we  see  that  E  £  H  (curl;  fi) 
satisfies 

(3.5)  /  ti-1^  x  E  -V  x£  -  (u)2e+iauj)E -ZdV  +  [  i  wT(nxE)-£dA  =  [  i  wJ-£dV 

Jn  Jr  is 

for  all  £  £  H (curl;  SI).  Existence  of  a  unique  solution  to  this  problem  follows  in  the  standard  way  [19]. 

Now  let  us  turn  to  the  inverse  problem.  Let  £  £  H( curl;  Si)  denote  a  solution  of  the  adjoint  Maxwell 
problem  (i.e.  satisfying  (2.7)).  Using  (3.5)  in  place  of  (2.6)  and  the  same  arguments  as  in  the  previous 
section  we  see  that  the  integral  equation  (2.8)  characterizing  J  is  modified  to  become 

(3.6)  iu>  J  J  ■  £  dA  =  J  nxE/i_1Vx£  +  icon  x  H  ■  £  dA. 

where  now  J  £  H-1/2( Div;  E)  and  £  €  H( curl;f2)  denotes  any  solution  of  the  adjoint  Maxwell  problem  (i.e. 
satisfying  (2.7)).  We  are  now  in  a  position  to  state  our  uniqueness  theorem  in  this  case. 

Theorem  3.1.  Suppose  J  £  f7_1/2(Div;  E)  and  E  and  H  satisfy  (3.1)-(3.f)  (i.e.  J  is  a  distribution 
supported  on  a  smooth  surface  E  in  Si).  Assume  also  that  E  is  a  given  surface  with  E  =  8B  and  B  is 
a  domain  in  ft  with  connected  complement.  Suppose  in  addition  that  the  electromagnetic  properties  of  the 
medium  in  B  are  constant.  Let  Eo  denote  the  support  of  J  on  E.  Then,  assuming  lo  is  not  a  Maxwell 
eigenvalue  for  the  adjoint  Maxwell  equation  in  B,  E0  and  J|s0  are  uniquely  determined  by  the  surface  data 
E  x  n  and  H  x  n  on  V. 

Remarks:  1)  This  theorem  rules  out  non-radiating  sources  on  smooth  surfaces  satisfying  the  conditions 
of  the  theorem.  We  note  that  for  the  Helmholtz  equation  Devaney  [8]  has  also  considered  the  problem  of 
non-radiating  surface  sources.  His  examples  show  that  for  more  general  surfaces  than  considered  here  (in 
particular  an  infinite  plane),  uniqueness  is  lost  since  non-radiating  sources  can  be  constructed. 

2)  If  a  >  0  in  B  the  eigenvalue  assumption  may  be  dropped. 

Proof.  By  virtue  of  the  assumptions  on  e  and  a,  the  electric  field  can  be  uniquely  continued  from  T  to 
Sl\B.  Thus,  from  the  point  of  view  of  uniqueness,  it  suffices  to  apply  (3.6)  on  dB.  Now  suppose  that  there 
are  two  currents  J\  and  J 2  in  H_1/2(Div;  E)  giving  rise  to  the  same  electromagnetic  fields  on  dB.  Then 

(3.7)  /  (Jr- J2)-£<M  =  0 

J  dB 

for  all  £  satisfying  the  adjoint  equation  (2.7)  (which  has  constant  coefficients  in  B). 

Let  £  be  a  Herglotz  wave  function  using  plane  waves  that  are  solutions  of  the  adjoint  Maxwell  system 
on  B.  It  is  shown  in  [5]  that  any  weak  solution  of  the  constant  coefficient  homogeneous  Maxwell  system 
can  be  approximated  in  the  H(curl;ft)  norm  by  a  Herglotz  wave  function  for  the  adjoint  system.  Given 
a  function  in  H-1/2(Curl;  E)  we  may  use  it  as  tangential  boundary  data  for  the  adjoint  Maxwell  problem. 
Using  the  fact  that  w  is  not  a  Maxwell  eigenvalue  for  the  adjoint  Maxwell  problem,  a  unique  solution  of  this 
boundary  value  problem  exists.  Approximating  the  resulting  solution  of  the  adjoint  problem  by  a  Herglotz 
wave  function  in  jH(curl;ft),  and  taking  the  tangential  trace  of  the  Herglotz  wave  function  shows  that  the 
tangential  trace  of  the  Herglotz  wave  functions  are  dense  in  H_1/2(Curl;  E)  (see  also  [19]).  Hence  (3.7)  shows 
that  Ji  —  J2  =  0  and  so  J  is  uniquely  determined  by  equation  (3.6).  We  have  thus  proved  the  theorem.  □ 
This  theorem  may  give  some  support  for  the  observation  that  rather  good  source  reconstructions  can  be 
obtained  from  MEG  data  assuming  the  source  is  on  the  brain  surface. 

4.  Uniqueness  for  dipole  sources.  The  final  case  we  shall  consider  is  where  the  current  source  is 
known  a  priori  to  be  the  superposition  of  a  finite  number  of  dipole  sources.  For  the  brain  these  dipoles  might 
represent  active  channels  in  the  nerve  membrane.  Thus  we  have  in  mind  that  a  finite  number  dipole  sources 
with  unknown  position  and  polarization  are  to  be  detected  in  the  interior  of  Si  (in  particular  the  sources  can 
not  lie  on  T).  This  is  not  a  restriction  in  practice  since,  when  detecting  sources  in  the  brain,  measurements 
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are  taken  on  the  skin  surface,  and  sources  are  located  in  the  brain  within  the  skull.  In  this  case  the  source 
takes  the  form 

i v 

J(x)  =  Yl  (*) 

j- 1 

for  some  N  <  oo  where  |Pj|  /  0  and  Zj  G  il  for  all  j  ( 5Z  is  the  delta  function  located  at  z).  We  assume  also 
that  the  source  points  Z\,  ■  ■  ■ ,  z^  are  distinct. 

In  order  to  establish  the  unique  determination  of  J  in  this  case,  we  need  an  additional  technical  restriction 
on  the  medium  near  the  source  points.  We  assume  that  each  source  point  is  located  in  a  region  of  constant 
electromagnetic  properties.  In  particular,  for  each  j,  there  is  a  ball  centered  at  Zj  such  that  e  and  p  are 
constant  in  the  ball.  This  assumption  is  used  first  to  justify  the  variational  formulation  of  the  inverse  problem 
(and  later  in  our  uniqueness  result) .  Since  a  weak  solution  of  the  Maxwell  system  is  a  strong  solution  away 
from  boundaries  and  discontinuities  of  the  coefficients,  we  know  that  the  solution  £  6  //(curl;  0)  of  (2.7) 
is  an  analytic  function  in  the  neighborhood  of  each  source  point  and  hence  (2.8)  still  holds.  In  the  case  of 
dipole  sources  (2.8)  gives 

N  _  f 

(4.1)  ^2/Pj  '  £(zj)  —  /  n  x  E  ■  p_1  V  x  £  +  iujn  x  H  ■  £ dA. 

i= i  Jv 

This  equation,  which  appears  already  in  [1],  will  form  the  basis  of  our  study  of  uniqueness  results  for  dipole 
sources  where  we  generalize  the  results  of  Ammari,  Bao  and  Flemming  [1]  to  more  than  one  dipole.  For  one 
dipole  in  free  space  this  theorem  was  also  proved  in  [13]  and  [3]. 

We  can  now  state  our  uniqueness  result. 

THEOREM  4.1.  Under  the  aforementioned  restrictions  on  the  electromagnetic  parameters  e  and  a  and  the 
source  points,  the  source  problem,  has  a  unique  solution  in  that  the  number,  position  and  non-zero  polarization 
of  the  point  sources  are  uniquely  determined. 

Proof.  The  first  part  of  the  proof  is  an  extension  of  the  proof  in  [1]  for  a  single  dipole.  We  start  by 
showing  that  the  number  and  position  of  the  dipoles  is  uniquely  determined.  Suppose  that  there  are  two 
sets  of  dipole  positions  S\  =  {z i,  •  •  • ,  z^}  and  S-2  =  {zi,  ■  ■  ■  ,zm }  giving  rise  to  the  same  measurements  on 
T.  By  assumption  the  points  in  Si  (respectively  52)  are  mutually  distinct.  If  Si  ^  S2  there  must  be  a  point 
in  one  set  not  in  the  other.  We  can  assume  z\  €  Si  and  Z\  S2.  By  unique  continuation,  the  electric  field 
can  be  continued  from  the  boundary  F  to  a  neighborhood  of  z\.  Let  Ei  denote  the  electric  field  continued 
to  the  domain  Q\S  1  and  E 2  denote  the  field  continued  to  the  domain  n\S2.  Then  E\  =  E2  in  a  punctured 
ball  around  z\. 

Now  suppose  that  x  approaches  z\.  The  field  E2  is  a  continuous  function  of  x  since  z\  is  assumed  to 
lie  in  a  ball  where  e  and  0  are  constant,  and  zx  is  not  a  source  point  for  E2.  Thus  \E2(x)\  — >  | jE?2 (^1 ) I 
as  x  — *  z\.  Considering  now  E\  there  are  two  possibilities:  either  pj  =  0  and  the  source  point  should  be 
dropped  from  Si  since  we  have  assume  the  the  polarization  at  each  source  point  is  non-zero,  or  Pj  ^  0.  In 
the  latter  case  |J57i  (cc) |  is  unbounded  as  x  — ►  zx  since  Z\  is  a  source  point  for  this  field.  This  contradicts  the 
previously  noted  equality  of  the  fields  and  so  Sx  C  S2.  Reversing  the  roles  of  Si  and  S2  shows  that  S\  —  S2. 

We  now  must  show  that  the  polarization  at  each  source  point  is  uniquely  determined.  We  know  already 
that  the  position  of  the  non-trivial  sources  is  uniquely  determined,  so  we  know  that  we  can  continue  the 
electric  field  to  a  punctured  ball  about  each  source  point.  In  particular  we  may  continue  the  field  to  the 
surface  of  a  ball  centered  at  a  source  point  containing  media  with  constant  electromagnetic  parameters  c 
and  cr.  Thus  we  are  faced  with  the  problem  of  knowing  n  x  E  and  n  x  H  on  the  surface  of  a  ball  containing 
a  single  source  point  at  it’s  center,  and  having  constant  electromagnetic  properties  inside.  If  B  denotes  the 
ball  and  dB  denotes  its  boundary,  the  variational  characterization  of  the  source  field  (see  (2.8)  becomes, 
after  moving  the  source  point  to  the  origin, 

(4.2)  p-m=  f  (nxil  p^Vxf-iwnxH  f)  dV 

JdB 

where  £  is  an  //(curl:  B )  solution  of  V  x  V  x  £  —  /c2£  =  0  in  B  and  k2  =  (w2ep  —  i crujp).  We  now  know 
that  £  =  r]exp(ikx  ■  d)  with  del3  such  that  |cf|  =  1  and  d  ■  77  =  0  is  a  solution  of  the  above  adjoint 
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Maxwell  system.  Using  this  solution  in  (4.2)  shows  that  p  f]  is  determined  for  any  17  and  hence  p  is  uniquely 
determined.  This  completes  the  proof.  Q 

5.  Measurements  on  a  portion  of  T.  In  real  applications  the  measurements  of  nx  E  and  nx  H  are 
only  made  on  a  subset  of  the  boundary.  We  now  assume  that  nx  E  and  nx  H  are  known  on  Fm  C  T  where 
rm  is  a  “measurement”  domain  that  may  be  multiply  connected.  In  this  case  a  variational  characterization 
of  J  can  still  be  derived.  Let  £  €  Hioc(curl;M3  \  Tm)  be  any  solution  of  the  adjoint  Maxwell  system  in 
K3  \  rm  together  with  the  adjoint  Silver-Miiller  radiation  condition  so  that,  in  a  weak  sense, 


(5.1) 

(5.2) 


Vx/i  1 V  x  E  —  (ui2e  —  iw0-)£  =  0  in  R3  \  F„ 

lim  ((V  x  £)  x  x  +  iwpoMO  =  0. 

|®J— >00 


We  then  have  the  following  lemma  giving  the  extension  of  (2.8)  to  this  case. 

Lemma  5.1.  Suppose  nx  E  and  nx  H  are  known  on  a  subdomain  Fm  of  F.  Then  provided  £  satisfies 
the  adjoint  problem  (5.1)-(5.2)  on  M3  \  Tm  then 


(5.3) 


iw  /  J  •  £dV  =  [  ( nxE ■  [p^V  x  £jT  +  i  umxH-  |£]T)  dA 
Jn  J  rm 


\  ft  in 


where  the  previously  introduced  tangential  jump  [-Jt  is  extended  to  a  smooth  function  on  ft  and 
the  obvious  way  so  that  for  x  €  T 

ImJt^)  =  lim  ((n(x)  x  u(x  —  hn( x)))  x  n(x)  -  (n(x)  x  u(x  +  hn(x)))  x  n(x))) 
h — ►() 

where  we  recall  that  n  is  the  unit  outward  normal  to  ft.  Analogues  of  equations  (3.6)  and  (4-1)  also  hold 
replacing  the  right  hand  side  of  the  respective  equations  by  the  above  right  hand  side. 

Proof.  The  proof  proceeds  similarly  to  the  derivation  of  (2.8).  Let  Bn  denote  a  ball  of  radius  R  with 
R  chosen  sufficiently  large  such  that  ft  C  Br.  Let  £  denote  a  smooth  test  function  on  Br  \  rm.  Then 
multiplying  (2.4)  by  £  and  integrating  by  parts  twice  we  arrive  at 

k o[j  £dV=[  E-hx  p-iV  x  £  -  (w2e  -  iwtr)£)  dV 
J  n  JBn\Tm  v  ' 

+  [  (n  x  E  ■  |/i-1  V  x  £jx  +  iom  x  H  ■  |£lr)  dA 
J  rm 

+  (nx  E  ■  p~lV  x  £  +  i  wn  x  H-l)  dA 
JdBR 

where  we  have  used  the  fact  that  J  =  0  outside  ft.  If  we  now  choose  £  to  satisfy  (5.1)  the  volume  integral 
term  vanishes.  The  integral  on  6Br  may  be  rewritten  as  follows: 

/  (nx  E  ■  /i_1V  x  £  +  ium  xH£)dA=  (e-  ((Pq1^  x  £)  x  n  +  iw^o£)) 

JdBR  JdBR  v  v  ' 


-iu(H  xn-E)-1)  dA. 


Using  the  fact  (see  [4])  that 

\E\ 


0(-j—r)  and  |(V  x  £)  x  n  +  iwpo^l  =  O 

\x\ 


X  ■ 


for  large  |x|  shows  that  as  R  — >  00  the  first  term  in  the  integral  on  the  right  hand  side  above  vanishes. 
Similar  estimates  show  that  the  second  term  vanishes.  This  proves  the  desired  equality.  □ 

It  is  now  possible  to  use  the  characterization  of  J  in  this  theorem  to  study  the  uniqueness  of  solutions  of 
the  inverse  source  problem  with  measurements  not  over  the  entire  surface  F.  In  addition  the  formula  could 
be  used  as  part  of  an  inversion  scheme. 

Obviously  volume  currents  cannot  be  reconstructed  using  (5.3).  However  the  proofs  of  uniqueness  in 
Sections  3  and  4  can  be  proved  provided  Tm  is  such  that  the  unique  continuation  property  holds. 
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6.  Conclusion.  Wc  have  examined  the  uniqueness  problem  for  the  inverse  source  problem  in  the  case 
where  the  background  medium  is  non-homogencous.  Wc  have  shown  that,  without  additional  a  priori 
information,  the  data  does  not  uniquely  determine  a  volume  source,  nor  does  it  even  determine  the  support 
of  a  volume  source  (these  results  are  well  known  for  a  constant  background) .  However  if  the  source  is  a  priori 
known  to  have  special  features,  in  particular  if  it  is  a  surface  current  on  a  known  surface,  or  a  a  current  due 
to  a  collection  of  dipoles,  we  show  that  the  data  does  uniquely  determine  the  current  density.  It  would  be 
useful  to  remove  the  many  technical  restrictions  regarding  the  electromagnetic  properties  assumed  in  this 
work. 

The  variational  characterization  of  the  current  distribution  in  each  case  can  also  be  used  as  part  of  inver¬ 
sion  scheme.  We  hope  to  test  this  approach  in  the  future,  in  particular  focusing  on  the  case  of  measurements 
on  a  portion  of  the  boundary. 
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